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On the MOND External Field Effect in the Solar System 
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Abstract In the framework of the MOdified Newto- 
nian Dynamics (MOND), the internal dynamics of a 
gravitating system s embedded in a larger one S is af- 
fected by the external background field E of S even if 
it is constant and uniform, thus implying a violation 
of the Strong Equivalence Principle: it is the so-called 
External Field Effect (EFE). In the case of the solar 



cessions will be independently estimated with different 
ephemerides it will be possible to repeat such a test. 
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motion through the Milky Way: it is orders of magni- 
tude smaller than the main Newtonian monopole terms 
for the planets. We address here the following ques- 
tions in a purely phenomcnological manner: are the 
Sun's planets affected by an EFE as large as 10 -10 m 
s~ 2 ? Can it be assumed that its effect is negligible for 
them because of its relatively small size? Does E induce 
vanishing net orbital effects because of its constancy 
over typical solar system's planetary orbital periods? 
It turns out that a constant and uniform acceleration, 
treated perturbatively, does induce non-vanishing long- 
period orbital effects on the longitude of the pericenter 
w of a test particle. In the case of the inner planets of 
the solar system and with E w 10~ 10 m s~ 2 , they are 
4 — 6 orders of magnitude larger than the present-day 
upper bounds on the non-standard perihelion preces- 
sions Atu recently obtained with by E.V. Pitjeva with 
the EPM ephemerides in the Solar System Barycen- 
tric frame. The upper limits on the components of E 



are E x < 1 x 10 



E v < 2 x 10- 



E z < 3 x 10" 14 m s~ 2 . This result is in agreement 
with the violation of the Strong Equivalence Principle 
by MOND. Our analysis also holds for any other ex- 
otic modification of the current laws of gravity yield- 
ing a constant and uniform extra-acceleration. If and 
when other corrections Ave to the usual perihelion prc- 
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The M Odified Newtonian Dynamics (MOND) (jMilgrom 
ll983allbU ch is a non-linear theory of gravity which under- 
goes departures from the standard Newtonian inverse- 
square law at a characteristic acceleration scale; it was 
proposed to explain certain features of the motion of or- 
dinary electromagnetically detectable matter in galax- 
ies and of galaxies in galactic clusters without resorting 
to exotic forms of still undetected dark matter. MOND, 
among other things, docs not fulfill the Strong Equiva- 
lence Principle in the sense that the internal dynamics 
of a gravitating system s of bodies does depend on the 
external background gravitational field E of a larger 
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Effect (EFE). 

Its effects have mainly been treated in a num- 

ber of ways in the context of galactic dynamics 

dSanders fc McGaughl20o3lFamaev et al.l2007aUAngus fc McGaugh| 
2008tlWu et al.ll2008l ). while they have received some at- 
tention in the case of t he Oort cloud in the solar system 
( Milgro nl ll983al . [r98l . [t is just the case to recall that 
it moves through the Milky Way at about 8.5 kpc from 
its center; its gravitational attraction at the solar sys- 
tem's location can be eval uated from the m agnitude of 
its cent rifugal acceleration (Milgrom|[l983b ) A ccn rj A , 



where ( Begcman et al.lll991 ) A = 1.27 x 10 
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m s 
Thus, 



is the MOND characteristic acceleration scale 
in principle, such an EFE should affect the inner dy 
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namics of the solar system. Concerning such an issue, 
a researcher active in MOND writes: "MOND breaks 
down the Strong Equivalence Principle. This means 
that the acceleration of solar system's bodies depends 
indeed on the background gravitational field and not 
only on the tidal field. As shown by Milgrom, even if 
the external field was constant (and the tidal force van- 
ishes), the internal acceleration would depend on the 
external field. Claiming that A con is irrelevant is only 
valid if the field equation were linear." He/she also 
adds that "for trans-Neptunian objects and planets, 
one can ignore the A ccn ." Another researcher working 
on MOND tries to go in deeper details by writing: "For 
the main planets, the acceleration is much larger than 
Aq (the order of magnitude of the EFE), and the effect 
is negligible [...] The EFE maintains a constant direc- 
tion in the planet revolution, and its effect cancels out. 
" Such statements are likely expressions of a widely 
diffuse belief about EFE in solar system. 

In this paper we deal with statements like those pre- 
viously quoted by showing, in a purely phenomenologi- 
cal way, that they are, in fact, not valid for the motion 
of the major bodies of the solar system referred to the 
Solar System Barycenter (SSB) frame usually adopted 
for planetary data reduction. In particular, we will ad- 
dress the following points 

• Does a constant and uniform acceleration E having 
a generic direction I and a magnitude of the order of 
E = A ccn « 10~ 10 m s~ 2 affect the motion of the 
major bodies of the solar system? 

• Are the effects of such an acceleration negligible for 
the Sun's planets, although they are, in principle, 
present? 

• Since E is constant and uniform, arc its effects can- 
celed out over planetary revolutions? 

We will demonstrate that 

• A constant and uniform acceleration does induce 
non-zero long-period, i.e. averaged over one orbital 
revolution, effects on the Keplerian orbital elements 
of a planet 

• By assuming E ss Aq, the resulting perihelion preces- 
sions of the inner planets are 4—6 orders of magnitude 
larger than the present-day limits on t he recently es - 
timated non-standard perihelion rates ( Pitjeval[2005 ). 



orbital revolution, can be modelled as 
E = E x i + E y j + E z k, 



(1) 



2 The impact of EFE on the dynamics of the 
inner planets 

Such an external field E, which would not be aligned 
with the internal Newtonian attraction TV, but it would 
be directed along a generic direction I fixed during one 



where E x ,E y ,E z are assumed constant and uniform, 
and i,j,k are the unit vectors of the usual SSB frame 
used to describe the planetary motions. The constraints 
we will obtain are, thus, of dynamical origin and model- 
independent. 

2.1 The perihelion precession induced by a constant 
and uniform perturbing acceleration 

From an observatio nal point of view, the astronomer 
E. V. IPitieval |2005l) has recently estimated, in a least- 
square sense, corrections Am to the standard Newto- 
nian/Einsteinian averaged precessions of the longitudes 
of the perihelia of the inner planets of the solar system, 
shown in Table [H by fitting almost one century of plan- 
etary observations of several kinds with the dynamical 
force models of the EPM ephemerides; since they do 
not include the action of any exotic modified model 
of gravity, such corrections are, in principle, suitable 
to constrain the unmodcllcd action of such a putative 
EFE accounting for its direct action on the inner plan- 
ets themselves. A preliminary, back-on-thc-envelope as- 
sessment of the impact of an additional acceleration E 
having the magnitude of Aq on the planetary perihelia 
can be done by simply taking the ratio of E to the av- 
erage speed of a planet along its orbit, given by v s=s na, 
where n — \J GM/a 3 is the un-perturbed Keplerian 
mean motion. For the Earth v « 3 X 10 4 m s _1 , so 
that w w E/na = 3 arcsec cty -1 , which is about 4 
orders of magnitude larger than the present-day uncer- 
tainty in the upper bound on any Earth's anomalous 
perihelion precession (Table H}. However, an accurate 
calculation must ultimately be done, also because our 
simple order-of-magnitude evaluation cannot tell us if, 
in fact, non- vanishing net secular perihelion precessions 
really occur under the action of a constant and uniform 
acceleration E. 



Table 1 

cty -1 (1 arcsec cty -1 = 1.5 x 1CP 15 s -1 ), to the standard 
Newton/Einstein perihelio n precessions of the inner plan- 
ets according to Table 3 of lPitieval (|2005h (Mercury, Earth, 
Mars). The result for Venus has been obtained by recently 
processing radiometric data from Magellan spacecraft (E.V. 
Pitjeva, private communication, 2008). The reference frame 
used is the usual Solar System Barycentric one. 



Mercury 


Venus 


Earth 


Mars 


-36 ± 50 


-4 ±5 


-2 ±4 


1±5 
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In order to precisely calculate the effects of eq. ([TJ 
on the orbit of a planet, let us project it onto the ra- 
dial (f), transverse (t) and normal (n) directions of the 
co-moving frame picked out by the three unit vectors 
|Montenbruck and GilllboOOl ) 



cos 51 cos u ~ cos / sin Q sin u 
r = | sin f2 cos u + cos / cos fl sin u 
sin / sin u 

- sin u cos f2 — cos / sin f2 cos u 

- sin Q sin u + cos / cos f2 cos u 
sin I cos u 

sin / sin Q 

- sin / cos f2 
cos/ 



(2) 



(4) 



where is the longitude of the ascending node which 
yields the position of the line of the nodes, i.e. the in- 
tersection of the orbital plane with the mean ecliptic at 
the epoch (J2000), with respect to the reference x axis 
pointing toward the Aries point T, / is the inclination 
of the orbital plane to the mean ecliptic and u = u> + / 
is the argument of latitude defined as the sum of the 
argument of the pericentre u>, which fixes the position 
of the pericentre with respect to the line of the nodes, 
and the true anomaly / reckoning the position of the 
test particle from the pericentre. Thus, it is possible to 
obtain the radial, transverse and normal components of 
the perturbing acceleration as 



E r = E • r, 



E,, 



E-t, 



E ■ h: 



(5) 



(0) 



(7) 



a straightforward calculation shows that they are lin- 
ear combinations of E x , E y , E z with coefficients propor- 
tional to harmonic functions whose arguments are lineal 
combinations of u, Q and /. E r , E t and E n must be in- 
serteH into trie riglit -liaurl-sirlp of the Gauss equations 
( Bertotti et al.ll2003l ) of the variations of the Keplerian 



orbital elements 



^eE r sin / + E t ( — 



dc 
dt 

dl 
dt 

dV. 

dl 

d,iu 

dl 



/I 



~ I 
na I 



E r sin f + E t cos / 



:(-:)]}■ 



na a/i"" 



- E n I — I cos u, 



na sin Iy/~1~ 



\/l - e 2 r / r\ I dQ 

— E r cos f -f- Ef I 1 H sin / - cos I , 

nae I \ p J J dt 



(3) d J± 

dt 



n E r 

na 



, / duj dVt \ 

^1 - e 2 h cos / 

V dt dt J 



(8) 

O) 
(10) 

(11) 
(12) 
(13) 



where e and A4 are the eccentricity and the mean 
anomaly of the orbit of the test particle, respectively, 
and p = a(l — e 2 ) is the semi-latus rectum. By evalu- 
ating them onto the un-perturbed Keplerian ellipse 



a(l-e 2 ) 
1 + e cos / 



(14) 



and averaging^ them over one orbital period Pb of the 
test particle by means of 



dt_ (1 - e 2 ) 3 / 2 
P b ~ 27r(l + ecos/) : 



;df, 



(15) 



it is possible to obtain the long-period effects induced 
byeq. ((2). 

In order to make contact. with the latest observa- 
tional determinations, we are interested in the averaged 
rate of the longitude of the periccntcr w = uj + f2; its 
variational equation is 



d-ca V'l" - 
dt nt 



-E T COS / + E t 1 H sin / + 2 sin 



(16) 



According to eq. (fT1)|) . it turns out that, after lengthy 
calculations, the averaged rate of w consists of a long- 
period signal given by a linear combination of E x .E y , E z 
with coefficients Cj, j = x,y, z of the form 



C 



F jk (e)cost;jk, j =x,y,z 



(17) 



where Fjk(e) are complicated functions of the eccen- 
tricity and £jfc are linear combinations of zu, O and /. 
In principle, they are time-varying according to 



VJ = TUg + ZOt, 

O = + f2*, 
/ = Io + it; 



(18) 
(19) 
(20) 



1 We used u = f because over one orbital revolution the pericentre 
ui can be assumed constant. 
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from a practical point of view, since their secular rates 
G7, 17, I are quite smaller, we can assume w ~ vjq, w 
£Iq,I ~ Iq in computing eos£fc, where zdo,QoiIo are 
their values at a given reference epoch (J2000); see on 



the WEB http://ssd.jpl.nasa.gov/txt/p_elem_tl.txt 



3 Constraints on EFE 



Such a constrain, which is a conservative one since we 
linearly added the uncertainties in E x , E y , E z to obtain 



SE < 



\EJ8EL 



\Ey\6Ey 



\EJ5Ez 



E 



(23) 



is four orders of magnitude smaller than A$. This sub- 
stantially confirms our preliminary, non-analytical eval- 
uation. 



Using the estimated extra-rates of Venus, Earth and 
Mars of Table [I] it is possible to write down a non- 
homogenous linear system of three equations in the 
three unknowns E x , E y , E z by equating the corrections 
Azu for the aforementioned planets to their predicted 
perihelion precessions due to E. It is 

Azj p — C px E x + CpyEy + C pz E z , p — Venus, Earth, Mars. 

(21) 

The values of the coefficients C P j are in Table fJJ 
they have been computed by using for vu, f2, / the val- 
ues at the reference epoch (J2000) by JPL, NASA 
( |http: / / ssd.jpl.nasa.gov / txt / p_elem_t 1 .txt [) , but it can 
be shown that Table [2] does not appreciably change 
if we vary them within the error bars of the EPM 2 . 
ephem erides retrievable, e.g., from Table 3 by iPitieva 
( 20081 ). As a result, we have for the Cartesian com- 



ponents of the perturbing acceleration E the figures 
quoted in Table [3] We have obtained them as follows. 
First, we have considered the eight different systems of 
the type of eq. (f_?Tj) obtained by taking the maximum 
and the minimum values of Aw for Venus, Earth and 
Mars according to their best estimates and uncertain- 
ties of Table [T] then, after solving them, we computed 
the means and the standard deviations of the three sets 
of eight solutions for E x , E y , E z . The upper bounds on 



E x . 

m s~ 2 , £ 2 <3X 10~ 14 
acceleration is 



,E y ,E z are E x < 1 x 10 

-,-2 



-15 



m s 



E y < 5 x 10~ 16 



m s , respectively. The total 



E= (0.6 ±3) x 10 



- 14 ms- 2 . 



(22) 



4 Interpretation of the results obtained 



How could we interpret the results obtained here? In 
fact, there is no contradiction with the fin dings in 



I Bekenstein k, Milgrom 1984 ; Milgrom 1986f) because 
they are valid in an inertial Galactocentric frame, for 
which the boundary condition at infinity yielding 



V<^ — > —E for r 



(24) 



where (p is the solution of the modified Poisson equation 
for the gravitational potential in the non-relativistic 
theor y of MOND proposed by (jBekenstein fc Milgrom 
198J), holds. See Figure [T] in which we have indicated 
the Newtonian monopole term with N and the Galactic 
attraction E with g cxt . 

It is not so in the SSB frame used for usual solar 
system's planetary data reduction because of the cen- 
trifugal acceleration which (almost) exactly cancels the 
Galactic gravitational attraction throughout SSB leav- 
ing just a small net tidal residual T, if any, whose mag- 
nitude can be easily evaluated as T < 10~ 14 m s -2 . See 
Figure O 

In other words, the outcome of the gravitational ex- 
periment represented by the planetary motion does de- 
pend on the velocities of the inertial frames in which 
it is studied, which are GC (at rest) and SSB (freely 
falling in the Milky Way). After all, this should not 
surprise too much, since MOND does, indeed, violate 
the Strong Equivalence Principle. 



2 In fact, EPM values of the planetary orbital elements the at 
(J2000) are not publicly available. 



Table 2 Computed values of the coefficients Cj, in s m~ , 
for Venus, Earth and Mars. For vj, Q, I entering cos£jfc 
the values at the reference epoch (J2000) have been used 
I http:/ /ssd.jpl. nasa.gov/txt/p_elem_tl.txtll. 





c 

^ x 


c 

V 


c z 


Venus 
Earth 
Mars 


0.0023 
-0.0005 
-0.0005 


-0.0052 
-0.0019 

9 x 10~ 6 


-0.0002 

5 x 10~ 10 
-1 x 10- 10 



Table 3 Cartesian components of the perturbing acceler- 
ation E, in m s~ 2 , from the means and standard deviations 
of the solutions of the eight different systems of equations of 
the type of eq. I|21|l obtained by combining the maximum 
and minimum values of Arh within their ranges according 
toTable[T](+ + +; + + -;H h; - + +; + --; h; - + 



E x (xlO- 15 ) 


Ey (XlO- 16 ) 


E z (xlO- 14 ) 


-0.3 ± 1 


2±5 


-0.6 ±3 
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Fig. 1 Sun-planet system s embedded in an external sys- 
tem S given by the Milky Way: Galactocentric frame. In 
it N (red arrow) is the Newtonian gravitational attraction 
of the planet by the Sun (the Newtonian gravitational at- 
traction of the Sun by the planet is not depicted), while 
E = g ext (light green arrow) is the gravitational attraction 
of the planet by the Galaxy whose magnitude amounts to 
about E ~ 10~ 10 m s~ 2 (in principle, it should also account 
for the gravitational attraction by other sources external to 
the Galaxy itself as the galaxy M31 Andromeda, etc.: how- 
ever, their acceleration is of the order of ~ O.OlAo). 



5 Conclusions 

In this paper we dealt with the problem of the Ex- 
ternal Field Effect in solar system in the framework 
of the MOdified Newtonian Dynamics. In particular, 
we analyzed the common belief that, in principle, also 
the internal motions of the Sun's planets are affected, 
in MOND, by a constant and uniform acceleration E 
having the same magnitude of the centrifugal accelera- 
tion of the solar system's revolution around the Galactic 
center A ccn « 10~ 10 m s~ 2 , and directed along a fixed 
direction I in space. It is believed that its magnitude is 
too small to affect the motions of the major bodies of 
the solar system, and that its net effects vanish because 
of its constancy during one planetary orbital revolu- 
tion. We found, instead, that such a kind of perturb- 
ing acceleration does, in fact, produce non- vanishing 
net long-period, i.e. averaged over one orbital revolu- 
tion, orbital effects on the motion of a test particle. 
Among them, of particular relevance are the perihelion 
precessions because it is possible to compare them to 
the corrections to the standard Newtonian/Einsteinian 
perihelion rates of the inner planets recently estimated 
with the EPM ephemerides in the SSB frame. As a re- 
sult, we obtained upper bounds on the components of 
E by finding E x < 1 X ICT 15 m s~ 2 , E y < 2 x 1CT 16 
m s~ 2 , E z <3x 1Q- 14 m s~ 2 , with E <3x ICT 14 m 
s~ 2 , i.e. 4 — 6 orders of magnitude smaller than A cen . 
This result can be interpreted in terms of a violation 
of the independence of the outcome of a gravitational 
experiment from the velocity of the frame in which it is 
performed, in according with the fact that MOND vio- 
lates the Strong Equivalence Principle. Our results are 
valid not only for EFE in MOND, but also for any other 
putative exotic modification of gravity yielding a con- 
stant and uniform net small extra-acceleration which, 
thus, can be no larger than 10~ 14 m s~ 2 in the solar 
system. It will be important to repeat such a test if and 
when other teams of astronomers will independently es- 
timate their own corrections to the standard perihelion 
precessions with different ephemerides. 
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Fig. 2 Sun-planet system s embedded in an external sys- 
tem S given by the Milky Way: heliocentric frame. In it N 
is the Newtonian gravitational attraction of the planet by 
the Sun (red arrow), while E is given by the gravitational 
tidal forces T (olive green arrow) exerted on the planet by 
the Galaxy whose magnitude amounts to about T < 10 -14 
m s -2 ; they account for both the uniform centrifugal ac- 
celerations due to the motion of Solar System through the 
Galaxy (light blue arrow) and the non-uniform gravitational 
attraction by the Galaxy throughout the extension of the 
Solar System (light green arrow). 
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